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ABSTRACT 

A  procedure  for  solving  the  plane  elastic  problem  for  edge  notches  in  a 
semi-infinite  region  is  presented.  A  modification  of  the  conventional 
Muskheli shvi li  technique  is  made  by  utilizing  a  mapping  function  which 
describes  both  the  physical  region  and  its  reflection  with  respect  to  the 
straight  edges.  An  effective  power  sdnes  development  in  the  parameter  plane 
is  then  described  and  illustrated  by  the  numerical  solution  for  semi -el lipti cal 
notches  in  a  semi -inf ini te  sheet  in  tension  at  infinity  parallel  to  the  edge. 
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LITERATURE  CITED 


INTRODUCTION 


Edge  notches  in  a  semi-infinite  region  identify  a  familiar  problem  class 
m  plane  elasticity*  The  difficulty  of  finding  a  representation  of  the  solu¬ 
tion  which  allows  a  systematic  consideration  of  boundary  conditions  on  both 
the  straight  line  and  notch  sections  of  the  geometry  is  a  challenging  matter. 
Several  particular  geometries  have  been  successfully  analyzed  by  a  variety  of 
rather  ingenious  methods.  Perhaps  the  most  familiar  is  the  original  Maunsell* 
solution  for  the  problem  of  a  semicircular  notch  in  a  semi - in f ini te  plate  in 
tension  parallel  to  the  edge.  Later,  Ling2  discovered  an  integral  represen¬ 
tation  for  this  problem  as  well  as  a  class  of  problems  corresponding  to 
circular  notches  or  mounds  and  straight  boundaries.  Edge  cracks  have  been  suc¬ 
cessfully  studied  by  several  investigators,  e.g.  ,  References  3  and  4.  Estimates 
of  the  maximum  stress  for  edge  notches  have  been  made  by  approximating  arguments 
by  Neuber. 5 


For  general  notch  shapes  a  more  uniform  approach  would  be  desirable.  One 
could,  for  example,  consider  the  mapping  technique  of  Muskhelishvili®  as  such 
an  approach.  However,  several  practical  difficulties  arise  in  the  application 
of  this  technique  to  this  problem  class.  Frequently  the  exact  mapping  function 
contains  branch  point  singularities  necessary  for  the  description  of  the  junc¬ 
tions  of  the  notch  and  the  straight  line  boundaries.  Effective  use  of  the 
exact  mapping  function  in  such  cases  is  limited,  as  the  problem  ceases  to  be 
one  of  linear  relationship.  The  alternative  of  polynomial  approximation  of  the 
exact  mapping  function  can  be  used,  but  frequently  a  suitably  accurate  poly¬ 
nomial  representation  is  awkward  to  find.  A  good  example  of  this  difficulty 
is  shown  in  the  recent  work  of  Mitchell^  in  his  solution  for  a  semicircular 
notch. 

In  this  paper  a  power  series  technique  is  presented  which  appears  to  pro¬ 
vide  an  effective  and  fairly  uniform  approach  to  single  edge  notches  in  a  semi¬ 
infinite  region.  Initially,  Muskhel i sh vi li ’ s  reflection  argument  is  used  to 
translate  the  problem  into  one  of  finding  a  suitable  analytic  stress  function 
defined  exterior  to  the  closed  region  bounded  by  the  notch  and  its  reflection 
with  respect  to  the  real  axis.  A  departure  from  the  conventional  mapping 
approach  is  then  made  by  introducing  a  mapping  function  defined  on  the  unit 
circle  and  its  exterior  to  describe  both  the  given  physical  region  and  its  re¬ 
flection  with  respect  to  the  real  axis.  A  power  series  representation  in  the 
parameter  plane  is  a  natural  consequence  of  this  formulation.  Furthermore,  in 
many  practical  problems,  the  explicit  occurrence  of  the  branch  point  singular¬ 
ities  in  the  mapping  function  referred  to  above  is  eliminated.  Finally,  the 
structure  of  the  solution  can  be  predicted  in  a  fairly  systematic  manner  by 
utilizing  most  of  the  consequences  of  the  unmodified  Muskheli sh vi li  theory. 


INITIAL  FORMULATION  OF  THE  PROBLEM 


It  will  be  assumed  that  the  mat 
S~  defined  in  the  z-plane,  Figure  1. 
assumed  to  be  an  arbitrary  continuou 
to  lie  on  the  real  axis,  y*0. 
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Figure  I.  SINGLE  EDGE  NOTCH  IN  A  SEMI-INFINITE  REGION 

In  the  initial  stage  of  the  formulation,  the  reflection  argument  of 
Muskhelishvili6  will  be  utilized.  The  essential  details  of  this  argument 
will  be  briefly  summarized  in  this  paragraph.  The  stress  components  in  carte¬ 
sian  coordinates  can  be  expressed  in  terms  of  two  analytic  functions  z)  and 
¥(z)  in  the  following  manner: 

+  vx  •  2  [*(*)  +  •(*)]  (1) 

+  2irxy  «  2[z$'(z)  +  ¥(z)],  (2) 

where  primes  denote  differentiation  and  bars  denote  complex  conjugates. 

In  the  present  problem,  the  region  of  definition  of  ®(z)  and  'Hz)  is  S”. 
On  the  other  hand,  using  the  extension  principle  of  Muskhel i shvil i ,  the  defi¬ 
nition  of  4>(z)  can  be  extended  into  S+  by  defining 

$(z)  *  -  t(z)  -  z$'(z)  -  flz)  for  z«S+,  (3) 

where 

f(z)  *  f(Tl.  (4) 

The  stress  function  Viz)  can  now  be  expressed  in  terms  of  the  extended  defini¬ 
tion  of  4>{  z )  as 

'f(z)  ■  -  $(z)  -4(z)  -  z$'(z)  for  zeS",  (5) 
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For  simplicity  of  presentation,  sel f- equi li bra  ting  load  systems  acting  on 
the  boundary  of  the  notch  will  be  considered  at  this  time.  Actually,  this  is 
a  minor  restriction  since  loading  at  infinity  or  along  the  straight  edges  can 
be  handled  directly  for  the  unnotched  half-plane  in  terms. of  elementary  func¬ 
tions  and  Cauchy  integrals,  e.g.,  Reference  6.  By  considering  the  superposition 
of  solutions,  it  is  clear  that  the  problem  difficulty  rests  in  solving  the 
cases  of  self-equilibrating  load  systems  acting  on  the  notch  boundary.  For 
self-equilibrating  load  systems, 

#(z)  and  ¥(z)  *  0(l/z2)  for  large  |z|.  (6) 

Thus,  from  (3)  it  is  evident  that  $(z)  is  analytic  in  both  S“  and  S+. 

The  character  of  $(z)  on  the  straight  edges  of  the  boundary  is  clear  from 

-  ir„v  =  *(z)  -  $(7)  +  (z-7)$'(z).  (7) 

If  the  straight  sections  along  the  real  axis  are  considered  as  load-free,  then, 
since  z-T’  *  0  on  the  real  axis, 

$“(x)  -  $+(x)  «  0,  (on  the  straight  boundary).  (8) 

The  notation  0~(x)  and  4>  +  (x)  denotes  the  value  of  $(z)  as  z—x  through  S“  and 
S+,  respectively.  From  (8)  and  the  previous  discussion,  4(z)  can  now  be 
considered  as  an  analytic  function  in  the  region  of  the  plane  exterior  to  the 
contour  C  =  C~  +  C+  where  C+  is  the  reflection  of  C-  with  respect  to  the  real 
axis.  Figure  1. 

Finally,  the  boundary  conditions  on  C"  must  be  considered.  For  the  pres¬ 
ent  purpose,  we  assume  that  the  normal  and  tangential  stress  components,  N 
and  T,  are  prescribed  on  C~.  Then,  if  a  denotes  the  angle  between  the  x-axis 
and  the  outward  normal, 

N  “  iT  *  $(z)  +  $(z)  -  e2la[(7-z)$'(  z)  -  $(z)  -  4>(z)]  for  zeC-.  (9) 

A  POWER  SERIES  FORMULATION  OF  THE  PROBLEM 

For  a.  power  series  representation,  it  is  obviously  desirable  to  deal  with 
a  parametric  region  consisting  of  a  circle  and  its  exterior  (or  interior). 
Therefore,  consider  an  auxiliary  complex  plane,  the  £-plane,  such  that  the 
unit  circle  £  =  cr  =  elS  and  its  exterior  are  mapped  into  C  and  its  exterior, 
respectively,  by  the  analytic  function 

z  -  *{£).  (10) 

A  departure  from  the  conventional  mapping  approach  is  made  here  since  both  the 
physical  region  and  its  reflected  image  are  described  by  the  auxiliary  plane. 

The  stress  functions  $(z)  and  f(z)  can  be  considered  as  analytic  functions  of 
the  parameter  £  since  z  and  £  are  related  by  the  analytic  function  (10).  The 
necessity  of  introducing  considerable  new  notation  can  be  avoided  by 
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designating  $  ( z )  =-$[w(  £ )  ]  as  $(£),  etc.,  which  leads  to  such  relationships  as 
'(z)=  *'(£)/«'(£),  etc.  Due  to  the  symmetry  of  C  with  respect  to  the  x-axis, 
*)=¥(£>  in  the  notation  described. 

The  boundary  condition  (9)  can  now  be  expressed  as  a  function  of  cr.  Due 
to  the  symmetry  of  C,  the  range  0<&<jr  will  be  considered  as  corresponding  to 
C+  and  the  range  Tr<9<2rr,  corresponding  to  C~.  The  applied  load  can  be  con¬ 
sidered  as  a  function  of  cr,  i.e., 

N  ~  iT  =  g(cr),  v<e<.2n.  (11) 


Thus, 


4>(cr)  +  cr)  -  e2l<x{  [co(cr)  -  co[a)]i  (cx)/a>'  ( cr )  -  $(cr)  -  $(cr)}  =  g(o-), 

TT<e<2n  (12) 


where 


e2ia  =  a2co‘(c r)/cu' (cr)'.  (13) 

On  the  unit  circle  itself,  it  is  possible  to  encounter  singularities  in 
$ ( £ )  either  by  load  distribution  or  as  a  result  of  geometric  irregularities 
in  the  notch  shape.  For  severe  types  of  singularities,  e.g. ,  poles  occurring 
at  crack  roots,  logarithmic  singularities  corresponding  to  point  loads,  etc., 
it  is  generally  possible  to  anticipate  and  introduce  the  proper  structure  as 
a  separate  part  of  the  representation.  The  remaining  arguments  will,  there¬ 
fore,  be  confined  to  the  cases  of  nonsevere  boundary  singularities,  i.e., 

$(£)  -  f  KOn,  I C ! >1  (14) 

n  =2  “ 

with  at  least  conditional  convergence  assumed  on  the  unit  circle  T. 


DETERMINATION  OF  THE  SERIES  COEFFICIENTS 

The  analyticity  of  $(£)  permits  the  expression  of  the  coefficients  in 
terms  of  the  well-known  integrals 


2*iAa  =  J‘*(l/£HC/£n+1,  n  *  2,  3,  ... 


(15) 


or 


2^ 


27rAn  ■  J  o"n$(cr) d£,  n  =  2,  3,  ••• 


(16) 


in  addition, 


277 


o  .  J  <r-n$(a)de,  n  -  1,  0,  -1,  -2 .  (17) 
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but  if  the  solution  exists  in  the  hypothesized  form,  the  conditions  (17)  ran 
be  considered  as  redundant. 


If  the  function  F(cr)  is  defined  as 

F{c)  =  fcia)  -  ♦(a)  =  0,  O<0<tt 

=  $(ct)  _  e 2ia  [  g(  cr)  -  <Mcr)  -  $(cr}]  -  4>(cr) 

+  [a>[cr)  -  co{cr)]  4'  [cr)/co'  (cr)  ,  v<0<2ir ,  (18) 

then  the  boundary  conditions  of  the  present  problem  are  equivalent  to 

F(o-)  =  0,  0<  S<  2w .  (19) 

The  Fourier  criteria  for  the  vanishing  of  F(cr)  is  thus 
Z.it  2w 


f  cr‘nF(cr)dfi  -  J 

cr‘nF(cr)d£  ■  0, 

n  =  0, 

±1,  ±2,  ... 

(20) 

0  v 

If  we  denote  by  Fj  , 

2n 

r>-l 

cr-iF(cr)d*, 

j  “  o, 

±1,  ±2,  ... 

(21) 

it  is  evident  that  computationally  (16)  and 

(17)  are 

equivalent  to 

Fj  *  0. 

j  -  2,  3, 

4 ,  ... 

(22) 

F.  -  0, 

j  -  1,  o, 

-1,  -2, 

•  •  • 

(23) 

respectively.  For  an  infinite  system,  the  redundancy  of  (17)  implies  the 
redundancy  of  (23). 


A  reasonable  plan  of  solution  conceptually  is  to  introduce  the  series 
(14)  into  the  set  of  conditions  (22).  This  leads  to  an  infinite  set  of  linear 
equations  in  the  unknowns  An,  which  can  be  solved  by  successive  truncations 
of  the  system.  Early  numerical  solutions  indicated,  however,  that  the  rate 
of  convergence  of  this  truncation  procedure  was  too  slow  to  be  practical.  The 
difficulty  shows  up  numerically  in  a  lack  of  dominance  of  the  principal  diag¬ 
onal  of  the  coefficient  matrix. 

Conceptual ly ,  the  difficulty  can  be  traced  to  the  series  truncation. 

Let  ♦ A( C )  be  defined  as  a  truncation  of  (14), 

M+l 

*.({)  -  Z  A  C * n .  (24) 

A  n*2  n 

In  general,  the  corresponding  function  F^(cr)  is  a  polynomial  (or  a  rational 
function)  of  a  involving  both  positive  and  negative  powers  of  cr  and  (19) 
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cannot  be  identically  satisfied  by  a  choice  of  the  M  available  coefficients 
An.  Furthermore,  the  Fourier  conditions  (23)  which  cannot  be  considered  as 
redundant  for  a  truncated  stress  function  are  ignored  in  the  plan  described 
above. 

A  combination  of  two  procedures  was  found  to  be  the  most  effective  com¬ 
putational  plan.  For  the  leading  coefficients,  the  system 

F,  =  0,  j  =  0,  ±1,  ±2,  ...,  ±P  (25) 

was  solved  for  several  values  of  P  with  M  =  2P  +  1.  A  rather  dramatic  improve¬ 
ment  in  convergence  of  the  leading  coefficients  was  found  by  this  plan.  In 
the  illustrative  problems,  the  first  twenty  coefficients  were  found  with  re¬ 
liable  accuracy  with  values  of  P<20.  On  the  other  hand,  a  critical  system 
size  is  reached  where  the  higher  ordered  coefficients  begin  a  rapid  growth  in 
inaccuracy  and  the  determinant  of  the  coefficient  matrix  becomes  very  small. 

This  can  be  anticipated  due  to  the  redundancy  in  the  limit  of  P  + 2  of  the 
equations  in  (25). 

The  procedure  used  to  calculate  the  higher  order  coefficients  consisted 
of  solving  the  set  of  conditions 

{F(o0}2  d6  -  0,  i  -  2,  3 . M  +  1.  (26) 

0 

This  set  of  conditions  corresponds  to  the  requirement  that  (19)  be  satisfied 
for  a  finite  set  of  coefficients  in  the  least  square  sense.  The  system  (26) 
eliminates  the  previous  stability  difficulties  for  larger  systems  and  can  be 
compactly  expressed  as  a  suitable  linear  combination  of  the  elements  Fj . 

Although  (26)  could  be  used  to  determine  the  entire  set  of  coefficients,  the 
leading  coefficients  converge  somewhat  more  slowly  than  in  (25),  and  the  max¬ 
imum  system  size  is  reduced  by  using  a  combination  of  the  two  procedures. 

SEMICIRCULAR  EDGE  NOTCH 

The  simplest  illustration  of  the  preceding  analysis  is  provided  by 
Maunsell's  problem,  i.e.,  a  semicircular  notch  in  a  semi-infinite  plate  in 
tension,  T,  parallel  to  the  edge,  Figure  2.  The  radius  of  the  semicircle  will 
be  chosen  as  the  unit  length,  and  polar  coordinates  (r,0)  are  introduced  in  the 
manner  shown.  For  illustrative  purposes,  the  mapping  rotation  will  be  retained, 
thus, 


The  stress 


functions 

4(C) 

YX(C) 


z  =  w( C )  =  C. 

will  be  denoted  by  4jU)  and  Y^(() 
=  T/4  +  4(C)  =T/4  ♦  T  2  At*" 

n  =  2  » 

-  -T/2  +  Y(C). 


where 


(27) 


(28) 
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Figure  2.  SEMICIRCULAR  EDGE  NOTCH 


Since  the  particular  solution  in  (28)  corresponds  to  cr x  =T  everywhere,  it  is 
clear  that  the  determination  of  $( £ )  and  y(i)  is  a  problem  type  with  the 
properties  assumed  in  the  previous  analysis.  Furthermore,  since  a  -  n  +  6, 

g(o-)  *  -T/2  d+cr2).  (29) 

From  symmetry  considerations,  the  series  coefficients  in  (23)  can  be  written 
in  the  form 

A2k  =  B2k>  k  =  1,  2,  3,  . . . 

A2k-1=  iC2k-l*  k  =  2,  3,  4,  ...  (30) 

where  and  02^.1  are  real. 

Corresponding  to  (18), 

F(cr)  =  0,  0<£<77 

=  $(cr)  +  4(cr)  -  (cr-cr- x)  ♦'  (cr)  +  a2[$(cr)  +  <J>(cr)] 

+  (T/2)  (1+ct2), 


tt<0<2jt. 


(31) 


Thus, 


Fo  =  wB2  -8  k?2  C2k- i/(  2k-3)  (2k+ 1)  +  rr/2 

Fi  =  0 

F2  =  82  C2k.  1/(2k-3)(2k-l)  +  w/2 
k  =  2 

F2j  =  2tt(  1-j  )  B2  j  +  (  2  j  -  1 )  7T  B2  j  _  2 

+  8  J2  C2k- i  (  2k-  1)  (  2j  -  1)/ [  (  2k-2j  )  2-  1]  [2k+  2j  -  3]  ,  j  =  2,  3,  ... 

00 

-iF2j.i  =  -(2j-l)C2j.rMj.l)^j.3-16  ^  B2k(2j-3)(j-l)/[(2j,3)2-4k2]  [2k-2j+l] 

-  4(j-l)/(2j-l)(2j-3),  j  =  2,  3.  ... 

F-2 j  =  7rB2j+2"8  J  C^jj.  1(2j  +  l)(2k-l)/[(2k+2j)2-l]  [2k-2j - 3] ,  j  =  1,  2,  ... 

‘iF-2j  +  r  7rC2j+l+16  B2kj(2j  +  l)/[(2k+2j)2-l][2k.2j-l]  ~4j/(4j2-l) , 

j-1,  2,  ...  / 


(32) 


The  conditions  corresponding  to  (26)  can  be  conveniently  expressed  in 
terms  of  Fn.  Since, 
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sr 


f”{F(o-) }zda  -  2  Y FU) 

rr  7T 


drier) 

Ba. 


de 


(33) 


we  find  from  the  even  and  odd  coefficients,  respectively, 

( 2  -  2  j  )  F_  2  j  +  (  2  j  +  1 )  F_  j  j  _  2  +  F2j_2  =  0,  j  =  1,  2 . M 

C2j  -l)Fi_2j  -  2j  F_1_2j  +  F2j.3  =  0,  j  =  2,  3,  ....  M  +  1.  (34) 

The  first  twenty  coefficients  were  found  to  five  decimal  place  accuracy 
by  using  (25)  for  several  values  of  P  <  20.  The  higher  order  coefficients 
were  then  calculated  by  solving  the  system  (34)  for  several  values  of  M. 

Forty  of  the  higher  order  coefficients  were  obtained  to  five  decimal  place 
accuracy  with  the  use  of  systems  in  which  M<50.  The  comparatively  rapid  decay 
of  the  coefficients  is  shown  in  Table  I. 


Due  to  some  recent  controversy  as  to  the  accuracy  of  Ling’s2  numerical 
results,  e.g.  ,  Reference  7,  the  stress  at  the  notch  root,  9  =  -tt/ 2,  was 
calculated.  In  terms  of  the  present  formulation, 
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Table  I.  STRESS  FUNCTION'  COEFFICIENTS  TOR 
SEMICIRCULAR  EDGE  NOTCH 


k 

B2k 

C2k  +  1 

k 

»2k 

C2k+1 

1 

•  0.79248 

-0.66283 

1 1 

-0.00  147 

-0.00074 

2 

♦  0.6  5206 

+0. 33908 

12 

-0.90110 

-0.09070 

3 

-0.04406 

+  0.05755 

13 

-0.00084 

-0. 00065 

4 

-0.0232S 

♦0.01765 

14 

-0.00065 

-0.  000.59 

5 

-0.01358 

+0.00616 

15 

-0.00051 

-0.00033 

6 

-C. 00850 

+0.00  199 

16 

-0.00040 

-0.90047 

7 

-0.00  561 

+0.00032 

17 

-0.00032 

-0.00042 

8 

-0.00385 

-0.00038 

18 

-0.0025 

-C. 009 38 

9 

-0.00272 

-0.00066 

19 

-0.00020 

-0.00034 

10 

-0.00198 

-0.00074 

20 

-0.000  16 

-0.000  30 

*  -*/2)  =  T  {1-4  [I  ( c2k  +  l-  B2k )  ( -  Dk]  }  *  (35) 

Using  the  coefficients  in  Table  I,  J  /T  =  3.0656  at  the  notch  root.  Allowing 
for  round-off  in  the  present  data,  it  would  be  reasonable  to  estimate  that 
3 X/T  -  3.065  ±  0.001  which  can  be  compared  with  Ling’s  result,  ~X/T  » 3. 065. 

Estimates  of  the  accuracy  of  the  solution  at  the  free  corner,  8  =0,  were 
made  from  the  data  of  Table  I  and  additional  higher  order  B„k  coefficients 
not  recorded  here.  Excellent  agreement  was  found  which  further  indicates  that 
reliable  peak  stresses  at  notch  roots  can  be  obtained  by  this  procedure,  par¬ 
ticularly  since  series  convergence  would  be  expected  to  be  more  rapid  at  such 
points. 


SEMI-ELLIPTIC  NOTCH  IN  A  SEMI-INFINITE  SHEET 

A  final  illustration  of  the  analysis  is  provided  by  a  semi-elliptic  notch 
in  a  semi -in fini te  sheet  in  tension  at  infinity  parallel  to  the  edge,  Figure 
3.  Only  an  approximation  of  this  problem  has  been  previously  available. 

Again,  the  straight  edges  of  the  sheet  are  assumed  to  lie  on  the  real  axis, 
y  *  0. 


An  appropriate  fc rm  of  the  mapping  function  is 

i  =  w(  O  =  C-BC*1,  0<B<1. 

The  ratio  of  the  axes  of  the  ellipse,  \  ■  a'b,  is  clearly 

K  ■=  (  1  -  B)  /  (  1  +  B) . 


(36) 

(  37) 


For  illustrative  purposes 
stress  functi on  $ ( £ )  and  ( ( ! , 


4 


the  analysis  was  carried  out  in  terms  of  the 
who  re 


6’U)  •  ♦  UVU),  *’(c)  ■  TUla'U!. 


(38) 
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The  preceding  arguments  carry  through  in  a  straightforward  manner  and  the 
boundary  condition  can  be  written  as 

<v cr)  $(cr)  a  u‘ (&)  $(cr)  +w(cr)4'(a)  -  u.'{cr)  <J>'  (a) 

-a;'(cr}  fco’icr)  g(cr)dcr,  *<9<2i t.  (39) 

-1 

If  $1^ £)  denotes  the  required  stress  function,  then 

♦  l(  O  =  (TAMt-Br1)  +  $(j)  (40) 

where 

4>U)  =  T  I  a.rn  +  i.  (41) 

n  =  2  n 

The  previous  argument  was  applied  to  the  function  cu'(  £)$(()•  Again  from 
symmetry, 

a2k  —  *  k  =  1 ,  2,  •  •  •  • 

a  2  ic  - 1  =  *c2k-l'  ^  =  2,  3,  ....  (42) 


10 


Since 


g(cr]  =  -(T/2)[l+  a2uj'  [a)/co,(cr)]  ,  (43) 

it  follows  that  in  the  boundary  condition  (39), 

f  co1  (  or)  g(  cr)  dcr  =  -  (  T/  2)  (  1  +  B)  (cr-a~  ‘ .  (44) 

Thu  s, 

f;  =  V  cr~  ]  {go1  (a)  <t>{  cr)  -  co'  {cr)  $  Icr)  +  $''  ( er'  [o;(  a)  -  co(  a )] 

J  v 

~co'  (a)  (  T/2)  (  1  +  B)  (cr-cr*  1 )  } d«9  ,  j  =  0,  ±  1 ,  ±2,  ...  (45) 

The  least  square  conditions  similar  to  (26)  can  then  be  written  as 
[2  +  B-2(  1  +  B)  j]  fU2j  +  [  2j  (  1  +  B)  -  1]  f_  2j  _  x 

*  ^  2  j  - 1  "  2  j  -  3  =  j  =  2,  ... 

[2(1  +  B)  j  -  (  1  +  2B) ]  f2_2j  -  [2(1+  B) j  -  (  2  +  3B) ]  f_2j 

-  T  2  j  -  2  ”^2j-4  =  j  =  2,  3,  ...  (46) 

The  numerical  analysis  was  then  carried  out  for  several  values  of  B  in 
the  manner  previously  outlined.  For  brevity,  only  the  peak  notch  stresses  for 
various  values  of  A.  »  a/b  wi  1 !  be  presented.  In  terms  of  the  present  formulation, 

^ax  "  °x(  £  =  *i}  =  T  U  +  4  |  [(  1  -  2k)b2k  +  2kczk  +  J  (-l)k/(l  -B)}.  (  47) 

Comparison  with  the  conventional  approximation  to  this  problem  was  made  by 
cal culating 

^  "  amax//'cr  max,  (48) 


where  cr*Bax  is  the  corresponding  peak  stress  for  an  elliptical  hole  in  an  infi¬ 
nite  sheet  with  a  corresponding  uniaxial  tension  at  infinity.  The  latter 
solution  is  well  known,  e. g. ,  Reference  8,  and 


=  T(  1  +  2\  - 1 ) . 


The  results  are  shown  in  Table  II  where  the  limiting  value  of  Q  for  A  =  0  is 
based  on  Koiter's  result.*  The  results  listed  can  be  considered  accurate  to 
within  0.1  of  one  percent.  In  the  calculation  of  the  coefficients  on  a 
digital  computer,  it  was  found  advisable  to  use  double  precision. 
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Table  II.  MAXIMUM  STRESSES  FOR  SEMI-ELLIPTIC  NOTCHES 


B 

\ 

cr  /T 

max' 

o*  *  /T 

max 

Q 

0.0 

1.0000 

3.065 

3.000 

1.022 

0.  1 

0.8182 

3.  540 

3.444 

1.028 

0.2 

0.6667 

4. 136 

4.000 

1.034 

0.  3 

0.5385 

4.910 

4.714 

1.042 

0.4 

0.4286 

5.948 

5.667 

1.050 

0.5 

0. 3333 

7.412 

7.000 

1.059 

0.6 

0.2500 

9.625 

9.  000 

1.069 

0.7 

C.  1765 

13.320 

12.333 

1.080 

0.8 

o.  mi 

By  interpolation  1.092 

0.9 

0.0526 

By  interpolation  1.106 

1.0 

0.0000 

Koi  ter1 s 

resul  t  “*  1.1215 
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